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Abstract 

A simple model for electron-vibron interactions on charged buck- 
minsterfullerene C6o n- > re = 1, ... 5, is solved both at weak and strong 
couplings. We consider a single H g vibrational multiplet interacting 
with t\ u electrons. At strong coupling the semiclassical dynamical 
Jahn- Teller theory is valid. The Jahn- Teller distortions are unimodal 
for 71=1,2,4,5 electrons, and bimodal for 3 electrons. The distortions 
are quantized as rigid body pseudo-rotators which are subject to ge- 
ometrical Berry phases. These impose ground state degeneracies and 
dramatically change zero point energies. Exact diagonalization shows 
that the semiclassical level degeneracies and ordering survive well into 
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the weak coupling regime. At weak coupling, we discover an enhance- 
ment factor of 5/2 for the pair binding energies over their classical val- 
ues. This has potentially important implications for superconductivity 
in fullerides, and demonstrates the shortcoming of Migdal-Eliashberg 
theory for molecular crystals. 

PACS: 31.30,33.10.Lb,71.38.+i,7420.-z,74.70.W 

1 Introduction 

The soccer-ball shaped molecule C 60 (buckminsterfullerene) and its various 
crystalline compounds have ignited enormous interest in the chemistry and 
physics community in past two years C 60 is a truncated icosahedron. 
From a physicist's standpoint, the charged molecule is fundamentally inter- 
esting, because the high molecular symmetry gives rise to degeneracies in 
both electronic and vibrational systems. Thus, the molecule is very sensi- 
tive to perturbations. In particular, electron-phonon and electron-electron 
interactions are expected to produce highly correlated ground states and 
excitations. 

Superconductivity has been discovered in alkali doped buckminsterful- 
lerenes A 3 C 60 (A=K,Cs,Rb), with relatively high transition temperatures 
(T c fa 20°-30°K). There are experimental indications that the pairing mech- 
anism originates in the electronic properties of a single molecule. The pair 
binding energy is a balance of electron-vibron interactions 0, |3|, f|] and 
electron-electron interactions ||. The relative contributions and signs of 
the two interactions is under some controversy. 

The electron-vibron school has identified certain five-fold degenerate H g 
(d-wave like) vibrational modes which couple strongly to the ti u Lowest Unoc- 
cupied Molecular Orbital (LUMO) 0, [|, |], §]. Varma, Zaanen and Raghava- 
chari || as well as Schluter et al and, more recently, Antropov et al proposed 
that these modes undergo a Jahn- Teller (JT) distortion and calculated the 
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induced pair binding energies at several fillings. They used the classical ap- 
proximation, and restricted their calculation to unimodal distortions (defined 
later). The general conclusion of this approach is that, while the calculated 
A is sizeable, one still requires a large reduction of the Coulomb pseudo- 
potential [i* in order to explain the highest transition temperatures. On the 
other hand, Gunnarsson et al independently estimate a large /x* ~ 0.4, i.e. 
there is no mechanism providing such a reduction. 

However, estimates of the electron-vibron coupling constant g do not 
justify the classical JT approximation. Cgo is estimated by frozen phonon 
calculations to be in the weak coupling regime g < 1 where quantum correc- 
tions are important. 

In this paper (Part I) we study the isolated Cg^ charged molecule. In par- 
ticular, we shall reconsider the same JT model, but diagonalize the quantum 
Hamiltonian for the full range of the coupling constant. We shall find that 
quantum corrections to the classical JT theory introduce novel qualitative 
features, and are quantitatively important for the pair binding energies. 

The quantum fluctuations involve interference effects due to geometri- 
cal Berry phases. Berry phases appear in a wide host of physical phenomena 
[0, B|. Here we find it in the context of a "Molecular Aharonov-Bohm (MAB) 
effect", originally discovered by Longuet-Higgins ||. The MAB effect has 
important consequences on the vibron spectrum. For example, it produces 
half-odd integer quantum numbers in the spectrum of triangular molecules 
[§, [7J, an effect recently confirmed spectroscopically in Na3 |[10|| . This kind 



of Berry phase is important also in scattering of hydrogen molecules [11 



Recently, it has been suggested that a geometrical Berry phase may be rele- 
vant in fullerene ions |12|, |13|, [L4j]. Here we show that Berry phases produce 
selection rules for the pseudo-rotational quantum numbers and kinematical 
restrictions which effect the pairing interaction between electrons. Although 
the semiclassical and Berry phase description is appropriate in strong cou- 
pling, the level ordering and degeneracies are found to survive for arbitrary 
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coupling, particularly in the weak coupling regime, which is closer to actual 
Ceo- For this reason we devote a large portion of this paper to the semiclas- 
sical theory, which helps to build physical intuition for further extensions of 
the model. 

This paper is organized as follows: In Section ^| the basic model is intro- 
duced. Section § calculates the JT distortions in the classical limit. Section f| 
derives the semiclassical quantization about the JT manifold. The geometri- 
cal Berry phases are calculated, and their effect on the semiclassical spectrum 
is obtained up to order g~ 2 . Section || describes the exact diagonalization 
results, and compares them to the semiclassical theory, and to weak coupling 
perturbation theory. The pair binding energies are determined in Section 
|6|. In Section |7| we summarize the paper and discuss our main result: that 
the effective pair binding energies are larger by a factor of 3 than the pair 
interaction energy in Migdal-Eliashberg's theory. In a following paper |15| 
we shall extend the model to all A g and H g modes with realistic physical 
parameters. This will allow us to explore the experimental consequences of 
the electron-vibron interactions. 



2 The Electron-Vibron Model 

The single electron LUMO states of C 60 are in a triplet of t\ u representation. 
We consider the H g (five dimensional) vibrational multiplet which couples 
to these electrons. t\ u and H g are the icosahedral group counterparts of the 
spherical harmonics {Yx m }l n __ 1 , and {^2A/}|/=-2 respectively. By replacing 
the truncated icosahedron (soccer ball) symmetry group by the spherical 
group, we ignore lattice corrugation effects. These are expected to be small 
since they do not lift the degeneracies of the L — 1, 2 representations. 
The Hamiltonian is thus defined as 

H =H° + H e ~ v , (1) 
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where, 

H° = hujJ2 (blM + 1) c Lcms ■ (2) 

M ' ms 

b M creates a vibron with azimuthal quantum number M, and $ ms creates an 
electron of spin s in an orbital Y lm . By setting |i-*6we can discard the 
second term. 

The H g vibration field is 

u(£l) = -^(Y 2 * M (Cl)bl + Y 2M {(l)b M ) , (3) 
where Cl is a unit vector on the sphere. The t\ u electron field is 

v>,(ft) = E Y lm (siy ms . (4) 

m=— 1 

The electron- vibron interaction is local and rotationally invariant. Its form is 
completely determined (up to an overall coupling constant g) by symmetry: 

H e ~ v oc g fdttu(£l)J2^l{&)M&) • (5) 

J s 

Using the relation 

Jdtt Y LM (Q)Y lmi (Q)Y lm2 (Q) oc (-1) M (L, -M\lm X ', lm 2 ) , (6) 



where (■ • •) is a Clebsch-Gordan coefficient fig] , yields the second quantized 
Hamiltonian 

H e ~ v = ^ghu £ OL + (-1) m &-m) 

s,M,m 

x (2, M\l, -m; 1, M + m)4 sCjtf+ros . (7) 

The coupling constant g is fixed by the convention of O'Brien, who studied 
first this kind of dynamical JT problem JT7I] . Representation (|?p is convenient 
for setting up an exact diagonalization program in the truncated Fock space. 



2.1 The Real Representation 

The semiclassical expansion is simpler to derive in the real coordinates rep- 
resentation. The vibron coordinates are 

6 2 



9m = -TfK E M ^ (4 + (-l) m M , (8) 



where 



M^ m ^Q = (2 sign(^)) 2 (<V, m + sign(/x)5 M _ m ) , 

M M)0 = ^o- (9) 
{q^} are coefficients of the real spherical functions 

/„(n) = ^M^n) 

V o m 

^i?e(F 2H (fi)) /i=l,2 

^r 20 (fi) a* = o . (io) 

^/m(y 2H (Q)) /i = -l,-2 
We also choose a real representation for the electrons 

4 _ J_(A 



C xs — \ C ls + c -ls 



V2 



^2 ^ " cLl ') 

4 = 4- (ii) 

Thus the Hamiltonian in the real representation is given by 
H = H° +H e ~ v 



% -V3q- 2 qo - V3q 2 ~V^q~i 



H e - V = <?^£(4,4>4) 



2 , 



■y/3qi VSq-i -2q 



( c 

Cys 

(12) 
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This form of the JT hamiltonian is well known [[17], 0] . Since the Hamilto- 
nian is rotationally invariant, its eigenvalues are invariant under simultaneous 
0(3) rotations of the electronic and vibronic representations. 



3 Jahn— Teller Distortions (Classical) 



In the classical limit, one can ignore the vibron derivative terms in (|i~2), and 
treat q = {q^} as frozen coordinates in H e ~ v . The coupling matrix in H e ~ v 
is diagonalized by [|18|1 : 



T-\w) 



z - v^r 
z + v^r 





T{w) 



(13) 



where 
T = 



/ cos ip smip 
— sin ip cos ip 

V o oi 



cos 9 sin 9 

1 
sin 9 cos 9 



( cos0 sin0 
— sin0 cos0 | . (14) 

V o 
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w = (<f),9,ip) are the three Euler angles of the 0(3) rotation matrix T. In 
the diagonal basis of (0), the electron energies depend only on two vibron 
coordinates: 





q(0) 



z 


Vo7 



(15) 



By rotating the vibron coordinates q to the diagonal basis using the L = 2 
rotation matrix , one obtains 



E 



(16) 



m,m' ,fi' 



where M^ m was defined in 
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By flT6|), and the unitarity of D and M, \q\ 2 is invariant under rotations 
of w. Thus, the adiabatic potential energy V depends only on r, z, and the 
occupation numbers of the electronic eigenstates n«, where J2i n t — n - 

V(z,r, [Ui]) = ^-{z 2 + r 2 ) + (ni(z - V3r) + n 2 {z + V3r) - n 3 2z s j . 

(17) 

V is minimized at the JT distortions (z n , f n , n^), at which the classical energy 
is given by 

E* =minV(z n ,f n ,ni). (18) 

The JT distortions at different fillings are given in Table I. We define 0, 9 
as the longitude and latitude with respect to the diagonal frame ( "principal 
axes") labelled (1,2,3) (3 is at the north pole), z, f parametrize the Jahn- 
Teller distortion in the real representation (|ToD , as 

(u JT (9,4>)) =^(3cos 2 £-l) +^sin 2 ^cos(20). (19) 

In Table I we present the values of the ground state JT distortions at all 
electron fillings. We see that electron fillings n — 1, 2, 4, 5 have unimodal 
distortions which are symmetric about the 3 axis, while n = 3 has a bimodal, 
about the 3 and 1 axes. The two types of distortions are portrayed in Fig. 
[I], we depict the distortions of ( |T9"D for the unimodal and bimodal cases. 

4 Semiclassical Quantization 

At finite coupling constant g, quantum fluctuations about the frozen JT 
distortion must be included. In order to carry out the semiclassical quan- 
tization, we define a natural set of five dimensional coordinates r, z, w. w 
parametrize the motion in the JT manifold (the valley in the "mexican hat" 
potential V) and r, z are transverse to the JT manifold, since V depends on 
them explicitly. The transformation q(r,z,w) was given in (|T6|), and was 
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derived explicitly in Ref. [|18j to be 

q 2 = z-y/3 sin 2 9 cos 20 + r— (1 + cos 2 9) cos 20 cos2-0 
2 2 

— r cos (9 sin 20 sin 2-0 

qi = z-v^sin 2$ cos0 — r- sin 2$ cos0 cos 2-0 

+ r sin # sin sin 2-0 
g = £-(3 cos 2 6* — 1) + f- sin 2 6* cos 2-0 

g_i = z-v^ sin 2$ sin — r- sin 26 1 sin cos 2-0 

— r sin 9 cos sin 2-0 

q_2 = z-V3 sin 2 6 1 sin 20 + r-(l + cos 2 0) sin 20 cos2-0 
2 2 

— r cos 9 cos 20 sin 2-0 . (20) 
The velocity in R 5 is given by 

q(r(t), z(t), w{t)) = d r qr + d z qz + <9 ro <f ■ w . (21) 

Using (^) and ([H]), we calculate the classical kinetic energy in terms of the 
JT coordinates. After some cumbersome, but straightforward, algebra the 
kinetic energy is obtained in the compact and instructive form: 



uj\ = — sin ip9 + cos if> sin 9<f>, 

u>2 = cos ip9 + sin ip sin 9<p, 

u 3 = if) + cos 9(f) 

(h,I 2 ,I 3 ) = ((V3z + rf,(V3z-r) 2 ,4r 2 ). (22) 

For finite JT distortions, we can identify the Euler angles terms as the kinetic 
energy of a rigid body rotator [19|], and the quantities Ii(z, f) as moments of 
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inertia in the principle axes frame 1,2,3. Thus, the Euler angles dynamics 
follows that of a rigid body rotator [ IB| . 

The unimodal and bimodal cases will be discussed separately. 



4.1 Unimodal Distortions 

For the unimodal cases (which we found for the ground states of n = 1, 2, 4, 5, 



r = on the JT manifold. The "moments of inertia" in ( 22|) are given by 
the tensor 

/l 0\ 

I = 3z 2 1 . (23) 
\0 0/ 

This corresponds to the rotational energy of a point particle on a sphere, 
which is described by the angles 8, <f>, and moment of inertia 3z 2 . Since axis 
3 has no "mass", its angular velocity is dominated by if). This implies that 
we must keep the term r 2 ip 2 but can discard the smaller mixed terms if?<j>. 
This yields 



The angular velocity if? couples to r 2 as in the kinetic energy of a three 
dimensional vector r parameterized by the cylindrical coordinates 

r = (r cos(2?/>), r sin(2?/>), z — z). (25) 

For |r| << z, the potential is simply 

V(r) ^|rf. (26) 

Thus, the semiclassical Hamiltonian of the unimodal distortion is 

jjuni ^ jjrot _|_ Jjho 
H rot = 

6z 2 

H ho = huj2(a\a 1 + l) , (27) 

7=1 A 
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where L is an angular momentum operator, and H are the three harmonic 
oscillator modes of r. The energies are given by 

E um = (^L(L + 1) + f> 7 + ~) j • (28) 
The rotational part of the eigenfunctions is 

= *Wfi) IKDfl, (29) 

where Cl = (6,<p) is a unit vector, and is the electronic adiabatic 

ground state. It is a Fock state in the principal axes basis. In terms of the 
stationary Fock basis |[n QS /]) where a = x,y,z, the adiabatic ground state is 

[n as ] 

Each overlap is a Slater determinant which is a sum of n products of spherical 
harmonics 

(WlWfl =E c MWn)W«)---Wn) - (si) 

M 

where C^j are constants. 

Now we discuss how boundary conditions determine the allowed values 
of L. A reflection on the JT manifold is given by 

Cl -> —Cl. (32) 

Spherical harmonics are known to transform under reflection as 

Ylw. ^( — l) L ^Lm • (33) 

Thus, by ( P0| ) and (0), the electronic part of the wave function transforms 

as 

IMa->(-i)1M-a- (34) 
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The reflection ( p2[ ) can be performed by moving on a continuous path on 
the sphere from any point to its opposite. (See Fig. |2[). It is easy to verify, 
using ( fT6| ) or (pOj) , that this path is a closed orbit of q e R 5 : 

01) -> g(-fi) = g(fi) . (35) 

Thus we find that the electronic wave function yields a Berry phase factor of 
(— l) n for rotations between opposite points on the sphere which correspond 
to closed orbits of q. In order to satisfy (p9"D) using the invariance of the 
left hand side under reflection, the pseudorotational Y^ m wavefunction must 
cancel the electronic Berry phase. This amounts to a selection rule on L: 

(-l) L + n = 1 . (36) 

Thus, the ground state for n — 1 and 5 electrons has pseudo- angular mo- 
mentum L = 1 and finite zero point energy due to the non trivial Berry 
phases. 

4.2 Bimodal Distortion 

The analysis of the bimodal distortions n = 3 proceeds along similar lines. 
The distortion obeys 

z = V3f . (37) 
From Eq.(|22"D we see that the kinetic energy is given by 



(38) 



where the inertia tensor is 

/4 

I = 2z 2 1 | . (39) 
\0 1 

The quantization of the pseudo-rotational part is the quantum symmetric 
top Hamiltonian. Fortunately, it is a well-known textbook problem (see e.g. 
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Ref. ^(J [16]]). The eigenf unctions of a rigid body rotator are the rotational 
matrices 

(40) 

where L,m,k are quantum numbers of the commuting operators L 2 ,L Z ,L X 
respectively. L z and L 1 are defined with respect to the fixed z axis and the 
co-rotating 1 axis respectively. The quantum numbers are in the ranges 

L = 0, 1, . . . oo 
m,k = —L, —L + 1, . . . L . (41) 

The remaining coordinates are two massive harmonic oscillators modes 

r = (r — f, z — z). (42) 

The semiclassical Hamiltonian is thus 

H bi pa H rot + H ho 

H = 4i^ L " 16i^ (L } ' 

H h ° = huj^{a\a 1 + \), (43) 

7=1 



and its eigenvalues are 

EU = ** (^L(L + 1) - ^~ 2 k 2 + E(n 7 + \) ) . (44) 

The rotational eigenf unctions are explicitly dependent on w as 

n o LM = D%l{w)X[\n ls ) VJ . (45) 

is 

4.3 Berry Phases of a Bimodal Distortion 

Unlike the unimodal case, in the bimodal case no single reflection fully clas- 
sifies the symmetry of the wavefunction. However, one can obtain definite 
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sign factors by transporting the electronic ground state in certain orbits. We 
define the rotations of it about principle axis U as Ci, which are schemati- 
cally depicted in Fig. |3|. The Berry phases associated with these rotations 
can be read directly from the rotation matrix T in Eq. fll4]) . For example: 
for tp—>ip + 7r (C 3 ), the states |1) and |2) get multiplied by (—1). 

Since transform as Yik under Ci, it is easy to determine the sign 

factors of the pseudorotational wavefunction. The results are given below: 

d : |1,0,2> 
C 2 : |l,0,2) w 
C 3 : |l,0,2) w 

q are coefficients in an L = 2 representation, and therefore are invariant 
under C\,C2,C 3 . Ci describe continuous closed orbits in R 5 . In order to 
satisfy ([K]) and using the degeneracy of E bl for k^—k, we find that 

L = odd , k = even. (47) 

In particular, the ground state of (|^) is given by L = 1, and k = 0. 

4.4 High-Spin Polarized Ground States 

It is possible to repeat the semiclassical analysis assuming that the spins are 
maximally polarized. These high-spin states are important, as they tend to 
prevail for strongly repulsive intra-level Hubbard U (Hund's rule) situations. 
In this case, we determine the JT distortions considering the Pauli exclusion 
between likewise spins. In Table II the JT distortions of the spin polarized 
ground states are listed. Our results for n=2,4 (S=l), and n=3 (5=3/2) 
cases are presented. The latter is trivial, since in that case n\ — n 2 — ns — 1, 
and therefore there is no JT effect at all. For n=2, (5"=1) there is unimodal 
distortion of z — —g which is smaller than the unpolarized ground state, 
and is equal to the distortion of the n = 5 case. Inspection of the orbital 



11,0,2), 
|l,0,2) t 
11,0,2), 



°1 • u m,k 
°2 • ^ m .k 



-l) k D. 



C, : D 



(L) 
m,k 



■l) L+k DtU 



(46) 
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energies t\ = €2 = —g , 63 = 2g provides a clear explanation for the identical 
distortions of the n=2 (5=1) and n=5 (S = |) cases, since in both cases e 3 
is occupied by a "spin up hole". 

Electronically, however, the two states are very different. First, we do not 
have a Berry phase for even number of electrons, as the individual contribu- 
tions from each of the two electrons cancel out. Second, there is a nonzero 
electronic orbital angular momentum. For example, the symmetry of the 
two-electron state prior to JT distortion is 3 P (i.e. 3 t\ u ), and so it remains 
following dynamical JT ||21|| . At finite coupling the two electrons in their 



ground state are still coupled in a 3 P electronic state, with L orb =l, where 
L or fy is the electronic orbital angular momentum, not to be confused with the 
pseudorotational quantum number L. Due to the absence of a Berry phase L 
must in fact be even, in contrast with the single electron case, and in agree- 
ment with Eq. (^). Thus although both cases have threefold degeneracies, 
they arise from different physical motion: purely electronic (for the n=2, 
S=l case) versus mixed electron-vibron motion (in the n=5, S=l/2 case). 

5 Exact Diagonalization 

The above semiclassical scheme gives a clear and intuitive picture of the 
behaviour of the system in a strong coupling limit. This limit is appropriate 
for describing, e.g., Na 3 [fLOR. However, in Ceo the actual range of the coupling 



parameter - g w 0.3 for a typical mode [p2| |15] - suggests that the electron- 
vibron coupling is actually in the weak to intermediate regime. 

Here we diagonalize the electron-vibron Hamiltonian (^) for single H g 
mode in a truncated Fock space. This approach yields accurate results unless 
the coupling strength is too large, and the higher excited vibrons admix 
strongly into the low lying states. We compare the results to the asymptotic 
large g expressions of the semiclassical approximation. The ground state 
energy for n = 1 has been previously computed in this fashion by O'Brien 
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H|. Here we present detailed results for all electron occupations, and also 
for the excitation spectra. 

Our basis is the finite dimensional Fock space of electrons and vibrons, 

l\n M ,n ms ) : N v < N max , £ n ms = n \ , (48) 

L ms ) 

where N v = J2m u m is the total vibron occupation. By gradually increasing 
N max , we have found empirically that accurate results can be obtained for g < 
N max /2, for levels with unperturbed energy below hu>N max /2. In particular, 
we have chosen N max = 5 (for n = 2,3) which yields an accuracy of better 
than 0.05huj for g < 0.6 and levels with N v < 1. The effect of truncation is a 
general upward shift of the levels, which gradually increases for higher excited 
levels. Level splittings and excitation energies are therefore less sensitive to 
the cutoff error. 

In Figures [|, || and | the ground state and a few of the excited states 
energies are plotted for one two and three electrons respectively. The four 
and five electron spectra are related to the two and one electron spectra by 
particle-hole symmetry. Energies are plotted as functions of g 2 . We compare 
the results to the semiclassical expressions (|28|) and (|44T ) for large coupling, 
and to second order perturbation theory at weak coupling. We discuss the 
different cases in detail, below. 

5.1 n = 1,5 electrons 

The ground state for one electron or hole in the t\ u shell is a threefold- 
degenerate state (all degeneracies given do not include spin) of the same 
symmetry: this fact is in complete analogy with what happens in the e <S> E 
coupled system, where the final dynamical JT coupled ground state has again 
E symmetry [^j]. Additional splitting of this ground state could occur via 



spin-orbit coupling, not included in the present treatment. Recent spectro- 
scopic data [E3] of Cg embedded in solid Ar confirm indirectly the presence 
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of the pseudorotational L = 1 ground state degeneracy, through direct obser- 
vation of spin-orbit splittings of about 30 and of 75 cm -1 for the t\ u ground 
state and for the t\ g excited electronic state m 1 eV above. The decrease 
of ground state energy is initially fast, and becomes gradually slower for 
increasing g. We shall return to this point in detail in 

As shown in Fig. [|, for large g, the n=l ground state energy correctly 
approaches the strong coupling limit 

1 9 3 1 
E ~ — g 2 + - + — , (49 

except for a small shift due, as mentioned above, to a finite-cutoff error. 
Above the ground state, there are families of excitations, corresponding to 
increasing values of N v . The lowest, for N v = 1, comprises 3 x 5 = 15 states, 
since for n—1, N v —1 there are just 3 electron states and 5 vibron states 
available. These states correspond to a direct product of a P (electronic) 
and a D (vibrational) manifold. As elementary angular momentum theory 
requires, they split into L = 3, 2 and 1 levels, which are found, in order of 
increasing energy. The splitting initially is proportional in g 2 , for small g, 
with significant deviations from linearity at g 2 m 0.2. As coupling increases, 
we note the slower downward trend of the even L states, than both the ground 
state and the associated "soft" odd-L excitations. This clearly reflects the 
Berry phase selection rule ( |36|) that no even L should appear among the low 
lying excited states in strong coupling. The lowest excitation from the ground 
state is L=l— > L—3, anticipating already at very weak coupling the strong 
coupling result that this excitation energy should fall fastest, and collapse as 
g^y. Unlike the L=3 state, the L=2 and L—l excited states do not show any 
tendency to collapse onto the ground state in the large g limit. Therefore 
they can be seen as modes involving essentially radial massive vibrations. 

The next group of excitations is for N v = 2, and comprise 3 x 15 = 45 
states. This multiplet splits into seven levels corresponding to L — 5, 3, 
1, 4, 2, 1 and 7. The lowest (L=5) level crosses two levels of the lower 
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(N v =l) multiplet in its downwards motion to become the second excited state 
above the L=3 level, eventually constituting the low energy odd-L rotational 
multiplet of the strong coupling picture. The same route is followed by the 
lowest level of N v = 3, which is an L=7 state. In fact, all the lowest split 
levels from each N v multiplet appear to have L—2N V + 1 and follow the same 
route. 

For N v = 2 we can similarly follow the movement with g of the L=4 level 
which decreases slowly towards the L=2 state from the lower N v — 1 to add 
to the group of massive radial vibrations. 

5.2 n = 2, 4 electrons 

Figure ||] has several features which contrast sharply with the one electron 
case. The N v = multiplet, has 15 two-electron states. The spin singlet sub- 
space constitutes of a 6-fold degenerate multiplet that splits into an orbital S 
and a D multiplet. As the semiclassical Eq. (^) suggests, the ground state 
and lowest excitations in the strong coupling limit have orbital degeneracies 
of even angular momenta. In fact, the lowest two among these states (L = 
0, 2) both come from the N v =0 multiplet, in contrast with the one electron 
case. The next pseudorotational level (L = 4) originates in the 6x5 = 30- 
fold degenerate N v — 1, spin singlet multiplet. Actually, at weak coupling it 
starts out being second in the ordering (L= 2, 4, 3, 2, 1, 0), but already at 
very small g it crosses the lower L=2 partner and approaches the pseudo- 
rotational asymptotic level. The convergence with increasing cutoff N max is 
worse than in the n=l case, which can be as due to larger JT distortions 
associated with two electrons. The spin triplet (5'=1) states of n=2 have not 
been plotted, as they behave in exactly the same fashion as the n=l states 
(see Fig. |]). This figure can be read in terms of n=2 S=l states simply by 
replacing the spin multiplicity label 2, as was in the case n—1, with 3. By 
comparison of Fig. ||] and Fig. § we notice that the low-spin 1 P state of N v =l 
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is exactly degenerate with the high-spin 3 D state in the same multiplet. This 
degeneracy seems accidental. 

5.3 n = 3 electrons 

For three electrons, the results are shown in Fig. |6|. The 8-fold degenerate 
N v = multiplet splits into two states characterized by degeneracies 3 and 

5 ( 2 P and 2 D). The ground state has the correct symmetry for an L—l, k=0 
state, which is predicted to be the ground state in the semiclassical limit. 
We also expect the lowest excitations to be classified as L=3, k—2 (14-fold 
degenerate), and L=3, k=0 (7-fold degenerate). In fact, three levels from 
the N v — 1 multiplet move down toward the ground state for increasing g. 
The one which moves lowest is 9-fold ( 2 G). In the g — > oo limit, it must 
therefore merge with the 5-fold levels from the N v =0 multiplet to produce 
the expected L=3, k—2 pseudo-rotator excitation. The next excitation of 
L=3, k=0 state can be identified as an asymptotic limit of the 2 F 7-fold 
degenerate state seen in Fig. |[ 

A remarkable feature of the n = 3 case is the presence in the N v — 1 
multiplet of a state (the 2 S) whose energy is independent of g\ This state is 
degenerate with the 5=3/2 state 4 D which has no JT distortion. 

6 Pair Binding Energies 

The pair energy for an average filling of n electrons is defined as 

U n = E n+1 + £ n -i - 2E n , (50) 

where E n are the fully relaxed ground state energies of n electrons. Formally, 
U is the real part of the two-electron vertex function at zero frequency. If 
this energy is negative for odd values of n, it means that electrons will have 
lower total energy if they separate into (n — 1) and (n + 1) occupations of 
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different molecules, rather than occupying n electrons on all molecules. For 
odd values of n, this is an effective pairing interaction often called "pair 
binding" in the literature ||. In Section |4] we found that for all odd n, the 
pair energies are negative, and given by the large g asymptotic expression 

£4=1,3,5 ~ -g 2 + i - ^ + (51) 

The first term is the classical energy. The second term is due to reduction 
of zero point energy along the JT manifold, since only radial modes remain 
hard. This term is independent of g and positive. The last term is due to the 
quantum pseudo-rotator Hamiltonian, and the Berry phases which impose a 
finite ground state energy associated with odd L for odd numbers of electrons. 
This term, although nominally small at large g, becomes important at weaker 



coupling. If ([51]) is extrapolated to the weak coupling regime the last term 
would dominate the pair binding energy. The exact diagonalization shown in 
Figure |7|, indeed shows a significant enhancement of the pair binding energy 
over the classical value in the weak coupling regime. 

In the weak coupling limit, we can obtain analytical expressions for U n (g) 
for g << 1 by second order perturbation theory. The unperturbed Hamil- 
tonian is the non interacting part H° with eigenstates (0). The perturbing 
hamiltonian is H e ~ v of Eq. (^), which connects Fock states differing by one 
vibron occupation. All diagonal matrix elements vanish, and the leading or- 
der corrections to any degenerate multiplet are of order g 2 . These are given 
by diagonalization of the matrix [p4fl , 

A nln> ms = ^n ms \H- EiQ) \ HQ H-\0,n' ms ), (52) 

in the degenerate 0-vibrons subspace. The sum implied by the inverse 
operator (E^ — H )^ 1 extends just to the N v =l states. The eigenvalues 
of yield the ground state energies and splittings for different electron 
fillings. These results, for all H g and also A g modes, and extended to the 
N v =l multiplet, will be discussed more extensively in |T5|. 
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Here we refer only to ground state energetics. In particular, using the 
perturbative expressions, we obtain, for a single H g , mode the small g pair 
binding energy 



n=l,3,5 



g 2 + o(g 4 



(53) 



The dependence strictly on powers of g 2 alone, with absence of all odd powers, 
is a consequence of the already mentioned AN V = ±1 selection rule of Eq. fl7|) . 
The origin of the 5/2 factor that characterizes the perturbative result ([53]) 
with respect to the classical pair binding energy (Table I) will be discussed 



in Ref. 15 



The molecular pair binding energy can be considered as an effective 
negative-U Hubbard interaction for the lattice problem, provided that the 
Fermi energy ep is not much larger than the JT frequency scale to. A 
mean field estimate of the transition temperature for the negative-f/ Hubbard 
model in the weak coupling regime is p5l l| 



ep exp 



(54) 



In Refs. || and |?B] , the results of Migdal-Eliashberg approximation for the 
superconducting transition temperature was given. Without the Coulomb 
pseudopotentials this approach yields 



V 



u;exp 



-6 9 



(55) 



By comparing ( |53D to flog) we find a striking discrepancy between the values 
of the effective pairing interaction: 



U =3V. 



(56) 



That is to say: in the weak coupling regime, the correct molecular calculation 
yields a pairing interaction which is three times larger than the results of 
Migdal-Eliashberg theory. 
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7 Discussion 



In this paper, we have solved the problem of a single H g vibron coupled to 
ti u electrons in a molecule. The model is too simplified for quantitative 
predictions for Cgo, but it contains interesting novel physics which will be 
important for further studies of this system. 

Semiclassically, a dynamical Jahn-Teller effect occurs. For n = 1,2,4,5, 
the molecule distorts unimodally, giving rise to a pseudo-angular momentum 
spectrum, plus three harmonic oscillators. For n — 3, there is a bimodal 
distortion, which generates a spectrum of a symmetric top rotator, plus two 
harmonic oscillators. The pseudo rotations are subject to non trivial Berry 
phase effects, which determine the pseudo-angular momenta L, and thus the 
degeneracies and level ordering of the low lying states. Strong Berry phase 
effects seem to survive even at moderate and weak coupling as shown by the 
exact diagonalization results. 

We find at weak coupling that the pair binding energy is a factor of 
5/2 larger than the classical JT effect, and a factor of three larger than the 
pairing interaction of Migdal-Eliashberg theory of superconductivity. This 
enhancement can be interpreted semiclassically as due to large zero point 
energy reduction of the pseudo-rotations. From the weak coupling point 
of view, this effect is due to degeneracies in both electronic and vibronic 
systems. 

Migdal's approximation neglects vertex corrections in the resummation 
of two-particle ladder diagrams. This is justified only in the retarded limit 
uj « ep. Here we have considered the opposite limit, where the molecular 
ground state energies are solved first, assuming that the JT relaxation time 
is of the same order, or faster than the inter molecular hopping time. In this 
regime, we have found therefore that Migdal's approximation substantially 
underestimates the pairing interaction, and T c , for these ideal molecular solids 



ZTfl . This large effect suggests that some of the enhancement is likely to carry 
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over to the real case of A3C60 metals, where electron hopping t and vibron 
frequencies are of similar strength. 

In Part II we shall consider a more realistic model which includes all 
important vibron modes of Cgo- We shall present quantitative predictions 
for the electron-vibron effects on the spectroscopy of C 60 ions. 
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Table I 



n 


S 




(711,712, n 3 ) 


E n /(hu) 


c„/(M 



1 
2 
3 
4 
5 
6 



1 

2 


1 

2 


1 

2 




(0,0) 
(0,0) 
(2s, 0) 

(-2 5 ,0) 

(-s,o) 

(0,0) 


(0,0,0) 
(0,0,1) 
(0,0,2) 
(1,0,2) 
(2,2,0) 
(2,2,1) 
(2,2,2) 


5 

_I„2 , 2 3 , J_ 
2# + 2 + 3<? 2 

-25 2 + | 

-lV + 1+3^ 

-2g> + | 

29 + 2 + 3<? 2 
5 
2 


V + i-* 

V + i-^ 



TABLE I. Semiclassical ground state distortions and energies for a single 
if s coupled mode of frequency to. n is the electron number, 5 is the total 
spin, z n , f n are the JT distortions, t!j is the occupation of orbital i, E n is the 
ground state energy and, and U n is the pair energy (Eq.(|5(|)). Energies are 
calculated for strong coupling to order g~ 2 . 



Table II 



n 


S 




(ni,n 2 ,n 3 ) 


E n /(hw) 


2 
3 
4 


1 

3 
2 
1 


(-5,0) 

(0,0) 

(5,0) (uni) 


(1,1,0) 

(1,1,1) 
(1,1,2) 


-19' + 1 

5 
2 

-V + - 

2 y ~ 2 



TABLE II. High spin ground state properties, in the same notation of 
Table I. 
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Figure 1: A polar representation of the Jahn- Teller distortions u JT (6,(f)), 
Eq. ([19|) . The distortion is measured relative to a sphere, (a) The unimodal 
distortion for the ground states of n = 1,2,4,5 electrons (b) The bimodal 
distortion for n = 3 electrons. 

Figure 2: Berry phase calculation for unimodal distortions. An path between 
reflected points on the unit sphere corresponds to a closed orbit in the five 
dimensional g-space. According to Eq. (|34]), such a path acquires a Berry 
phase of (—1)™ from the n-electron wavefunction. 

Figure 3: Berry phases calculation for the bimodal distortion (n = 3). w 
are the three Euler angles which rotate the principal axes of the bimodal 
distortion. denote rotations by 7r around corresponding axes. On the 
right we depict the electronic Berry phases associated with the three closed 
orbits in g-space, given by Eq. (f46[). 

Figure 4: Exact spectrum for one electron as a function of the square 
electron- vibron coupling constant g 2 . The vibron occupations are truncated 
at N max = 5. The semiclassical energies (Eq. (p8|) ) are drawn by dashed lines 
for the lowest three pseudorotational multiplets (n 7 =0, L—l, 3, 5). The unit 
of energy is the vibron quantum hou. 

Figure 5: Exact spectrum for two electrons (S = 0). The semiclassical 
energies (Eq. fl2"Bp) are drawn by dashed lines for the lowest three pseudoro- 
tational multiplets (n 7 =0, L=0, 2, 4). The two-electron S = 1 spectrum is 
the same as for n = 1, S = |. 

Figure 6: Exact spectrum for three electrons (S = |). The semiclassical ener- 
gies, Eq.([P]), are drawn by dashed lines for the lowest three pseudorotational 
multiplets, n 7 =0, (L, &) = (!, 0), (3,2) (3,0). 
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Figure 7: Pair binding energy U (thick solid line), compared to weak coupling 
perturbation theory for g « 1 (dotted line) and semiclassical theory for 
g » 1 (dashed lines). U n is found to be the same for n — 1,3,5 electrons. 
The Migdal-Eliashberg approximation V (thin solid line) is also drawn for 
comparison. g C(i0 is the range of physical coupling strength for Ceo- 
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